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' Abstract 

^ : 

T— I . The black-hole binary Cygnus X-1 was observed for 17 ks with the Suzaku X-ray observatory in 2005 

CO ' October, while it was in a low/hard state with a 0.7-300 keV luminosity of 4.6 x 10'^^ erg s~^. The XIS and 

\ HXD spectra, spanning 0.7-400 keV, were reproduced successfully incorporating a cool accretion disk and 

^-H ' a hot Comptonizing corona. The corona is characterized by an electron temperature of ~ 100 keV, and two 

\ optical depths of ~ 0.4 and 1.5 which account for the harder and softer continua, respectively. The disk has 

OO . the innermost temperature of ~ 0.2 keV, and is though to protrude half way into the corona. The disk not 

' only provides seed photons to the Compton cloud, but also produces a soft spectral excess, a mild reflection 

^ , hump, and a weakly broadened iron line. A comparison with the Suzaku data on GRO J1655— 40 reveals 

\^ ' several interesting spectral differences, which can mostly be attributed to inclination effects assuming that 

\ the disk has a flat geometry while the corona is grossly spherical. An intensity-sorted spectroscopy indicates 

■ that the continuum becomes less Comptonized when the source flares up on times scales of 1-200 s, while 
the underlying disk remains unchanged. 
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1. Introduction 

Luminous soft X-ray radiation of accreting stellar-mass 
black holes (BHs) has generally been explained as thermal 
emission from optically-thick (in particular "standard") 
accretion disks (Shakura & Sunyaev 1973; Makishima 
et al. 1986; Dotani et al. 1997; Remillard & McChntock 
2006), which are expected to form around them under 
rather high accretion rates. In contrast, their hard X-ray 
production process is much less understood, even though 
intense hard X-ray emission characterizes black-hoc bina- 
ries (BHBs) among various types of compact X-ray sources 
in the Milky Way and Magellanic clouds. 



Indeed. BHBs often emit a major fraction of their ra- 
diative luminosity in the hard X-ray band, in the form 
of spectral hard-tail component if they are in so-called 
high/soft state, or as the entire power-law (hereafter PL) 
like continua if they are in so-called low/hard state (here- 
after LHS) which appears under relatively low accretion 
rates. Furthermore, the hard X-ray emission (partic- 
ularly in the LHS) involves another interesting aspect, 
namely the long-known aperiodic variation over a wide 
frequency range (e.g., Oda et al. 1971; Oda 1977; Nolan 
et al. 1981; Miyamoto et al. 1991; Pottschniidt et al. 
2003; Remillard & McChntock 2006). These spectral and 
timing studies are not limited to stellar-mass BHs, since a 
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fair fraction of active galactic nuclei (AGNs) are also con- 
sidered to be in a state analogous to the LHS of BHBs, 
emitting luminous hard X-rays (e.g., Mushotzky, Done & 
Pounds 1993). 

In hard X-rays, a BHB in the LHS emits a roughly 
PL shaped continuum with a photon index F ~ 1.7 (e.g., 
Remillard & McChntock 2006), but the spectrum grad- 
ually steepens toward higher energies. This has been ex- 
plained in terms of unsaturated Comptonization of some 
soft seed photons, in a cloud of hot Maxwellian electrons, 
or a "corona" , with a temperature of Tg = 50 — 100 keV 
and an optical depth of r ~ 1 (Shapiro et al. 1976; Sunyaev 
& Triimper 1979; Sunyaev & Titarchuk 1980; Gierhnski 
et al. 1997; Di Salvo et al. 2001; Frontera et al. 2001a; 
Zdziarski & Gierhnski 2004; Ibragimov et al. 2005). As 
a result, the Comtpon y-parameter, y = ArkTo/mcC^ , also 
take a value of ~ 1, where k is the Boltzmann constant. 
Too is the electron mass and c is the light velocity. The 
electrons will be cooled by transferring their energies to 
the photons, while heated via Coulombic collisions by ions 
that acquire nearly a free-fall temperature (e.g., Zdziarski 
et al. 1999). The Comptonizing cloud may be identi- 
fied with radiatively inefficient inner accretion flow (e.g., 
Ichimaru 1977; Narayan & Yi 1995). 

Superposed on the dominant continuum characteriz- 
ing the LHS, we often observe a spectral hump at en- 
ergies of 20-40 keV due to "reflection" of the primary 
X-rays by some cool materials (e.g., Lightman & White 
1988; Inoue 1989), a fluorescent Fe-K emission line (Basko 
1978; Kitamoto et al. 1990), and an apparently "smeared" 
Fe-K edge (Ebisawa et al. 1994). These features can be 
explained by assuming that a fraction of the primary hard 
X-rays are reprocessed by cool materials located near the 
BH, presumably an optically-thick part of the accretion 
disk (Done et al. 1992; MarshaU et al. 1993; Ebisawa et al. 
1996; Gierhnski et al. 1997). This part of the disk may at 
the same time produce an X-ray excess which we often ob- 
serve at the softest end of the spectrum (e.g., Balucinska 
& Hasinger 1991), and furthermore, may provide the 
Comptonizing corona with the seed photons (e.g., Shapiro 
et al. 1976). 

In spite of general success of this currently popular in- 
terpretation of the LHS in terms of Comptonizing hot 
clouds and cool materials, details of the scenario still 
need to be elucidated. It is yet to be confirmed that the 
seed photons are provided by a cool accretion disk (e.g., 
Zdziarski et al. 1999), rather than by thermal cyclotron 
emission in the Compton corona itself (e.g.. Kantach et al. 
2001). The geometry and size of the Comptonizing cloud 
is not well known, except for loose constraints, includ- 
ing; that the cloud should be large enough for the com- 
pactness parameter to remain within a reasonable range 
(Lightman & Zdziarski 1987), but should not be too large 
to be able to produce the hard X-ray variation on a ~ 1 
s time scale, and that it should have a sufficient verti- 
cal height so as to efficiently illuminate the reprocessing 
matter. Still unsettled is the relation between the hot 
corona and the putative cool disk; whether the cool disk 
is truncated at a large radius (Zdziarski et al. 1998; Done 



& Zycki 1999; Zdziarski et al. 2004), or it penetrates deep 
into the hot corona (Young et al. 2001; Miller et al. 2006b). 
Yet another fundamental issue is how to understand the 
violent hard X-ray variability in the context of thermal 
Comptonization scenario (e.g., Miyamoto et al. 1987). 

To address these issues associated with the LHS of 
accreting BHs, a number of spectral studies have been 
conducted so far, by fitting observed broad-band spec- 
tra with theoretical spectral models (e.g., Gierhnski et al. 
1997; Zdziarski et al. 1998; Dove et al. 1998; Frontera 
et al. 2001a). Similarly, extensive studies in the frequency 
domain have been conducted to characterize the random 
variations (e.g., Miyamoto and Kitamoto 1989; Negoro 
et al. 1994; Nowak et al. 1999b; Gilfanov et al. 2000; 
Revnivtsev et al. 2000; Remillard & McClintock 2006). 
However, the former approach is limited in that differ- 
ent modelings often degenerate. The latter approach has 
also a limitation, in that it may not constrain details of 
the Comptonization scenario without invoking particular 
models (Kazanas et al. 1997; Hua et al. 1997). 

We analyze broad-band (0.7-400 keV) Suzaku data of 
Cygnus X-1 (Cyg X-1; Oda 1977; Lian and Nolan 1984) 
acquired in 2005 October. With the wide-band sensitiv- 
ity, Suzaku is expected to provide novel diagnostics of this 
BHB, the first celestial object suspected to be a black 
hole (Oda et al. 1971). To overcome the above difficul- 
ties, we compare the Suzaku spectra of Cyg X-1 with 
those of GRO 1655—40, searching the two spectra for dif- 
ferences attributable to their different system inclination, 
i - 45° of Cyg X-1 (Abubekerov et al. 2004) and i ~ 70° 
of GRO 1655-40 (Green et al. 2001). We hence refer to 
the Suzaku results on GRO 1655—40 by Takahashi et al. 
(2007a) as Paper I. We also analyze spectral changes asso- 
ciated with the random variations on a time scale of ^ 1 s 
(section 4) , to see how the Compton parameters are chang- 
ing; this will be complementary to the study of long-term 
spectral changes (e.g., Ibragimov et al. 2005). We employ 
a distance of 2.5 kpc (Margon et al. 1973; Bregman et al. 
1973) to Cyg X-1. The errors refer to statistical 90% con- 
fidence limits. 

2. Observation and Data Processing 

2.1. Observation 

The present Suzaku observation of Cyg X-1 was per- 
formed on 2005 October 5, from UT 04:34 through 15:11 
for a net exposure of 17.4 ks, as a part of the ini- 
tial performance verification program of Suzaku (Mitsuda 
et al. 2007). Suzaku carries 4 sets of X-ray telescopes 
(Serlemitsos et al. 2007), each coupled to a focal-plane X- 
ray CCD camera called XIS (X-ray Imaging Spectrometer; 
Koyama et al. 2007) operating in the energy range of 0.2- 
12 keV. Three (XISO, XIS2, and XIS3) of the four XIS sen- 
sors use front-illuminated (FI) CCDs, while XISl utilizes 
a back-illuminated (BI) one thus achieving an improved 
soft X-ray response. The Hard X-ray Detector (HXD; 
Takahashi et al. 2007c; Kokubun et al. 2007) covers the 
10-70 keV energy band with Si PIN photo-diodes (HXD- 
PIN), and the 50-600 keV range with GSO scintiUation 
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counters (HXD-GSO). 

Because the optical axes of the XIS and the HXD arc 
shghtly (~ 3'. 5) offset, the observation was carried out 
with the source placed at the center of the XIS field 
of view. The HXD was operated in the nominal mode 
throughout the observation. Because the source was very 
bright, the XIS was operated with the "1/8 window op- 
tion" , in which a smaller field of view of 17'. 8 x 2'. 2 is read 
out every 1 s. The editing mode of the FI CCDs was set 
to 2 X 2, while that of the BI CCD (XISl) was 3 x 3. In 
the present observation, however, XISl with the BI CCD 
chip suffered severe telemetry saturation, for more than 
60% of the total exposure. 

During the observation, the averaged source count rates 
(excluding backgrounds) of XISO, HXD-PIN, and HXD- 
GSO were ~ 300, ~47, and ~20 c s~^, in the energy bands 
of 0.5-10 keV, 10-70 keV, and 70-400 keV, respectively. 
However, the true XISO counts would be higher, if the 
data were free from telemetry overflows and event pile up 
(subsection 2.2). 

2.2. XIS Data Selection 

We retrieved the XIS data which were prepared via ver- 
sion 1.2 processing. The XIS data were further screened 
based on the following standard criteria: the XIS GRADE 
should be 0, 2, 3, 4, or 6; the time interval after an exit 
from the South Atlantic Anomaly should be longer than 
500 s; and the object should be at least 5° and 20° above 
the dark and sunlit Earth rim, respectively. 

In order to conduct spectral and timing studies, we uti- 
lize XIS events extracted from a region which is defined 
to fit in the "1/8 window". Since the Suzaku attitude 
is known to wobble on a time scale of several tens min- 
utes, the region needs to be a little smaller than the win- 
dow size. Considering the magnitude of the wobbling, we 
defined the region as an intersection between a rectan- 
gle of 70 X 1024 pixels and a circle with a radius of 3'. 3 
(200 pixels), both centered always on the image peak. For 
this purpose, we calculated the image center every 200 s 
and re-defined the event extraction region also every 200 
s. Thus, the location of the region on the CCD chip is 
time dependent. The background was entirely negligible 
throughout the observation (less than 0.1% of the signal 
at 8 keV considering the Galactic ridge emission), so was 
not subtracted from the data. 

Because Cyg X-1 was so bright that we need to care- 
fully remove the effects of telemetry saturation and photon 
pile up. When the number of detected events exceeds the 
telemetry allocation to the XIS, the excess events are dis- 
carded without being transferred to telemetry. This phe- 
nomenon is called telemetry saturation. Since the event 
buffer is cleared every 8 s, and the XIS exposure was 1 s 
in the present observation, we excluded the telemetry- 
saturated data in the following way. For each 8-s win- 
dow, we calculated an average count rate over the first 4 
s which are usually free from the telemetry saturation. If 
a subsequent exposure in the same 8-s window showed a 
count rate lower than this average by more than 2 stan- 
dard deviations, we judged that the telemetry was satu- 



rated, and discarded this and the subsequent exposures 
in the same window. This typically occurred in the last 
one or two exposures in the 8-s window. Because these 
telemetry-saturated time bins arc regarded as dead times, 
the source count rate remains approximately intact, al- 
though the exposure is reduced. 

After filtering out the telemetry-saturated data, we fur- 
ther removed the effects of photon pile up. Although the 
Suzaku X-ray telescopes have relatively wide point spread 
functions, photon pile up becomes noticeable at the im- 
age center when observing a point source brighter than 
^ 100 cts per exposure. The average count rate of Cyg 
X-1 during the present observation was ~3 times higher 
than this limit. We therefore accumulated XIS events 
by excluding circular regions of various radii at the im- 
age center, and examined the obtained series of spectra 
for shape changes. As a result, we have decided to ex- 
clude an image-center region with a radius of I'.O (60 
pixels). Although this means that the event extraction 
region has a complicated shape (intersection of an annu- 
lus and a rectangle), we use standard ancillary response 
files (ae_xi0_xisnom4_20060615 . arf , etc.) defined for a 
circular region of 4'.0 radius. This is because the point 
spread function of the telescope has little energy depen- 
dence, and hence the spectrum would not change except in 
normalization. We confirmed this inference using the data 
of a moderately bright point source, SS Cyg observed on 
2005 November 2, and found that the spectral normaliza- 
tion should be multiplied by a factor of ~ 0.092 compared 
with that of PIN, when excluding the central 1' of the 
image; this is taken into account in the subsequent spec- 
tral analysis. We also use the standard response matrices 
(ae_xi0_20060213.rmf , etc.) for the XIS. 

2.3. HXD Data Selection and Background Subtraction 

The HXD data utilized here were prepared via ver- 
sion 1.2 processing. We further screened the HXD events 
by imposing conditions that more than 500 s should be 
elapsed from an exit out of the South Atlantic Anomaly, 
that the target should be > 5° above the Earth rim, and 
that the geomagnetic cutoff rigidity should be higher than 
8 GV. 

Since the HXD has neither imaging capability nor an 
offset counter to measure instantaneous non X-ray back- 
ground (NXB), we need to estimate the NXB and subtract 
it from the raw data. Accordingly, we employed the PIN 
background modeling by Watanabe et al. (2007), known 
as "PINUDLC (bgd_a)" model ^, which combines actual 
HXD-PIN data acquired while pointing to night Earth 
(Kokubun et al. 2007). Thus, we produced a large num- 
ber (10.0 times as large as the actual) of fake non X-ray 
background events, which would just be detected during 
the present observation if we had an imaginary NXB mon- 
itor detector. By processing this fake event files exactly 
in the same manner as the actual data, we can estimate 
NXB contributions in the PIN data, and subtract them 
from PIN light curves as well as from PIN spectra, with a 

^ http:/ /www. astro. isas.jaxa.jp/suzaku/analysis/hxd/hxdnxb/ 
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typical accuracy of 5% or better. The cosmic X-ray back- 
ground entering through the PIN field of view is entirely 
negligible in the present study. 

Subtraction of the HXD-GSO background employs a 
different method, developed by Fukazawa et al. (2007) and 
called "LCFIT (bgd_d)" model . In this case, the whole 
GSO energy band is divided into 32 energy bins, and a 
light curve of each energy bin is produced over a long 
period. Because each light curve has characteristic repe- 
tition patterns with periods of one orbital revolution, one 
day (due to the relation with respect to the South Atlantic 
Anomaly), and 51 days (due to orbital precession), it can 
be expressed as a function of time using a relatively small 
number of adjustable parameters. The parameters can 
be determined by analyzing these light curves over a long 
period, typically one month. We can thus derive the ex- 
pected NXB counts in individual energy bins at a given 
time during the present observation, and subtract them 
either from light curves or spectra. The GSO data are 
hence analyzed using the same energy binnings as used in 
the NXB modeling. When integrated over ~ 1 day, this 
method reproduces the GSO background within a sys- 
tematic error of ~ 2% (Takahashi et al. 2007b). Since 
the PIN and GSO background models both utilize PIN 
upper-discriminator hit rates, their basic time resolution 
is typically 4 s or longer, depending on the telemetry rate. 
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Fig. 1. Light curves of Cyg X-1 binned into 64 s, obtained in 
the present Suzaku observation performed in 2005 October. 
The ordinates are in units of ct s^^. From top to bottom, 
panels refer to those obtained with XIS2 (0.7—2 keV), also 
with XIS2 (2-10 keV), with HXD-PIN (10-60 keV), and with 
HXD-GSO (60-200 keV). The HXD backgrounds were sub- 
tracted. The periodic data gaps are due to spacecraft passages 
through the South Atlantic Anomaly and the Earth occulta- 
tion. 

2.4- Light Curves 

Figure 1 shows XIS2, HXD-PIN, and HXD-GSO light 
curves of Cyg X-1 from the present observation. The HXD 



backgrounds (PIN and GSO) were subtracted in the way 
as described in section 2.3, and the results were corrected 
for dead times. In contrast, the XIS background, which is 
completely negligible, is inclusive. The count rates indi- 
cated by these XIS light curves are by an order of mang- 
nitude lower than those actually observed, because of the 
image-center exclusion which was needed to avoid pile-up 
events (subsection 2.2). On time scales of a few minutes 
to a day, the source thus varied randomly typically by 
~ ±10%. 

2.5. Time- Averaged Spectra 

We show in figure 2 the time-averaged Suzaku spectra of 
Cyg X-1. Like in section 2.4, we corrected the HXD spec- 
tra for the dead times, and subtracted the backgrounds 
using the method described in section 2.3. In contrast, 
the XIS2 spectrum is here and hereafter not corrected for 
the normalization reduction (down to 0.092 of the original 
level), due to the exclusion of the central 1' of the images 
(subsection 2.2). For comparison, we superpose the mod- 
eled background spectra of HXD-PIN and HXD-GSO. 

As seen in figure 2, we have significantly detected Cyg 
X-1 nearly across the full energy band of Suzaku. In par- 
ticular, the source signals in HXD-PIN exceed the NXB 
by more than a factor of ~ 5, even at the highest energy 
of ^ 70 keV. Therefore, the systematic background uncer- 
tainty of HXD-PIN is completely negligible compared to 
photon counting statistics. Even in the GSO range, the 
source signal is so strong that it can be detectable up to 
^ 400 keV, above which the systematic background errors 
start dominating. 




Energy (keV) 

Fig. 2. Time-averaged XIS2 (black), HXD-PIN (red), and 
HXD-GSO (green) spectra of Cyg X-1, shown after subtract- 
ing the background but not removing the instrumental re- 
sponses. The XIS2 spectrum is accumulated after removing 
the telemetry overflow and and cutting off the image center 
(see text). The modeled non X-ray background of HXD-PIN 
is indicated in magenta. That of HXD-GSO is shown in cyan, 
together with its 5% and 3% levels. 



http: / / www.astro.isas.jaxa.jp/suzaku/analysis /hxd / gsonxb / 
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3. Analysis of the Time Averaged Spectra 

This section is devoted to the study of time-averaged 
Suzaku speetra of Cyg X-1. Because of tlie strong teleme- 
try saturation (section 2.1), we do not use the data from 
XISl (BI). Of the three FI CCD cameras comprising the 
XIS, only the data from XIS2 are employed in our spectral 
analyses, while those from XISO and XIS3 are utilized in 
section 4 for a different purpose; even if we incorporated 
data from the other XIS sensors, systematic errors would 
dominate before data statistics are improved, because Cyg 
X-1 is so bright for the XIS. 

3.1. Fitting with conventional models 

In order to roughly characterize the time-averaged 
broad-band spectra derived in subsection 2.5, we simul- 
taneously fitted simple empirical models to the XIS2, 
HXD-PIN, and HXD-GSO spectra, so that the results 
may be compared with other observations of Cyg X-1, as 
well as those of other similar objects. Here and hereafter, 
the overall model normalization is allowed to differ be- 
tween the XIS and HXD (mainly to absorb uncertainties 
introduced by the removal of telemetry saturated data 
and the image-center exclusion), but constrained to be 
the same between HXD-PIN and HXD-GSO. As the 
detector responses, we use ae_xi2_20060213 . rmf 
and ae_xi2_xisnom4_20060615.arf for XIS2, 
aeJixd_pinxiiioin_20060814.rsp for PIN, 

while aejixd_gsoxinom_20060321 . rsp and 

aeJixd_gsoxiiiom_20070424.arf for GSO. The last 
GSO file is the correction factor described in Paper I, 
introduced so as to make the Crab Nebula spectrum 
described by a single PL model with F = 2.1 ^. We 
discard the XIS data in the 1.7-1.9 keV range to avoid 
complexity of the Si-K edge modeling, and those in the 
> 8 keV band which is possibly affected still by the photon 
pile-up effects. The HXD-PIN data below 12 keV are also 
discarded in the fitting, to avoid response uncertainties 
there (Kokubun et al. 2007). We utiHze the HXD-GSO 
data in the 70-400 keV energy band, because of the poor 
signal significance (figure 2) and relatively large response 
uncertainties, in > 400 keV and < 70 keV, respectively. 
We assign a systematic error of 1% to all the spectral 
bins, in order to absorb uncertainties in the instrumental 
calibration, including in particular those associated with 
the detector responses and background estimation. The 
effect of contamination on the optical blocking filters of 
the XIS is modeled as additional absorption (Koyama 
et al. 2007). We also trimmed the energy offset by ~ 
17 eV, considering the current uncertainty with the XIS 
energy scale in the 2x2 mode. 

As the simplest attempt, we first fitted an absorbed sin- 
gle PL model simultaneously to the three sets of spectra. 
The fit indicated F^ 1.77 as the best estimate, and an ab- 
sorbing hydrogen column density of A^h = 3.5 x 10^^ cm~^. 
However, the fit was far from acceptable, with /v = 14.3 
for V = 360 degrees of freedom. As seen in figure 3b, the fit 



leaves several noticeable discrepancies, including a promi- 
nent data deficit above ~ 100 keV, and a data excess in 
0.9-1.5 keV. Obviously, the former is due to the high- 
energy cutoff, while the latter is due to the presence of 
soft excess, both of which have been established through 
previous observations. 



10 

1 

0.1 

0.01 

2.0 




.£ 0.7 



1.0 



I 

(3 0.7 

1.0 

0.7 



—P 1 1 1 — P I I I I 1 1 1 P — I I I I I P 1— 

( a) d i 5ltbb-i-cLitoff pl+retlection+gaus 5 
XIS X PIN 



■ \ 1^ 




(b) power-law (r=].77) 



H\ 1 1 I I I I ll| 1 1 I I I I ll| 1 r-p 

(c) diskbb+cLitotfpl 



I I I I ll| 1 — h 



(d) di s kb b+c uto rip I +re tlec tion+ga us s 
I ■ I I I > ' t > i i I I I I t t t it\ I 



12 i 10 20 50 100 2(M) 400 
Energy (kcV) 

Fig. 3. Simultaneous fittings to tlie time-averaged XIS2, 
HXD-PIN, and HXD-GSO spectra of Cyg X-1, us- 
ing various combinations of conventional models. (a) 
Background-subtracted and response-inclusive spectra (the 
same as in figure 2), compared with the best- fit 
model consisting of a disk-blackbody, a cutoff pow- 
er-law, a Gaussian, and a reflection component; i.e., 
wabs* (diskbb+cutof f pl+pexrav+gau) . Individual model 
components are separately drawn, (b) Data-to-model ratios 
for an absorbed single-PL fit. (c) The same as panel b, but 
a disk-blackbody component is added, and the power-law is 
replaced by a cutoff power-law model, (d) The data-to-model 
ratios corresponding to the fit in panel a, namely, obtained 
by further adding to the model of panel c a Gaussian and a 
reflection component. 

To express the high-energy turn over, we next re- 
placed the PL component with a cutoff power-law model 
(cuoffpl; a power-law times an exponential factor). To 
account for the soft excess, we also added a disk-blackbody 
component (diskbb: Mitsuda et al. 1986; Makishima et al. 
1986) to our fitting model. The diskbb model is a simple 
non-relativisitic approximation to the X-ray spectrum in- 
tegrated over a standard accretion disk, and is considered 
to provide an approximately (within ~ 10%; Watarai et 
al. 2000) correct innermost disk radius when an appropri- 
ate correction factor is incorporated (Kubota et al. 1998). 
Thus, the model is wabs* (diskbb+cutof f pi) in the xspec 
terminology; here, we utilize version 11 of the xspec spec- 
tral fitting package ^, together with the photoelectric ab- 
sorption model, wabs, by Morrison & McCammon (1983). 



http: / /www. astro. isas.jaxa.jp/suzaku/analysis/hxd/gsoarf/ 



ttp: / /www. heasarc.gsfc.gov/docs/xanadu/ 
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Then, the fit was hiiproved drastically to jv = 3.7 for 
V ~ 357, though not yet acceptable. The exponential cut- 
off energy was obtained as Tcut '^-^ 170 keV, and the nominal 
PL photon index decreased to F ^ 1.45. The innermost 
disk temperature of the diskbb component was obtained 
as Tin ~ 0.5 keV, and its innermost radius as r^n ~ 13 
km at a distance of 2.5 kpc (section 1), without any cor- 
rection for the color hardening factor, or inner boundary 
condition, or inclination (see Makishima et al. 2000). We 
regard these diskbb parameters still unphysical. 

In figure 3c, the data to model ratios from the 
wabs* (diskbb+cutof f pi) fit show a slight broad hump 
over 20-40 keV (though not very clear in this presen- 
tation), indicative of reflection from a cool thick ma- 
terial. In addition, there is clear evidence of Fe-K 
line at about 6.4 keV. Accordingly, we added a re- 
flection component (pexrav; Magdziarz & Zdziarski 
1995) and a Gaussian to our model, thus constructing 
a model wabs* (diskbb+cutof fpl+pexrav+gau) , and re- 
peated the fitting. In using pexrav, the reflector was as- 
sumed to have an inclination i = 45°, and its incident 
photons to have the same spectral shape and normaliza- 
tion as those of the cutof fpl continuum (rel_ref 1< 
technically). The Gaussian was allowed to have a free 
centroid and a free width. Then, the fit further improved 
to /v = 2.4 for V — 353, implying that the reflection 
is signiflcantly present as reported by many previous ob- 
servations. The model gave the cuoffpl parameters as 
F 1.74 and Tcut > 1000 keV, and the diskbb parame- 
ters as Tin ~ 0.3 keV and nn ^ 45 km. The reflector was 
found to have a moderate solid angle as ri/27r ~ 0.5, and 
the Gaussian centroid and its width were obtained as 6.43 
keV and a < 0.2 keV respectively. In figure 3a, we show 
this fit where contributions of the individual model com- 
ponents are separately drawn. The data-to-model ratios 
are given in figure 3d. Thus, the reflection signal is in- 
ferred to carry ~ 30% of the continuum at 30 keV. 

In this way, the empirical model flttings have roughly 
quantified the time- averaged Suzaku spectra of Cyg X-1, 
and extracted their basic features. The fit implies a 0.7- 
300 keV flux (after absorption by A^h) of 5.4 x 10~* erg 
cm~^ s~^. Assuming an isotropic emission, and employ- 
ing the distance of 2.5 kpc (section 1), the absorption- 
removed 0.7-300 keV luminosity becomes 4.6 x 10'^^ erg 
s~^, which is typical of this object. Although the fit still 
remains unacceptable, and hence errors associated with 
the model parameters cannot be evaluated, further im- 
proving the model would be of little meaning, because 
cutof fpl is a mere empirical approximation to the spec- 
tra expected from unsaturated Comptonization processes. 
Instead, we proceed to more physical modelings, using 
theoretical Comptonization codes. 

3.2. Fitting with Comptonization models 

Among various Comptonization codes currently avail- 
able, here we employ that of Poutanen & Svensson (1996), 
called compPS in xspec, because of the following three ad- 
vantages: it correctly takes into account the relativistic 
effects in the Compton scattering kinematics; it allows us 



to use a diskbb model as the seed photon source; and 
it has a built-in function to allow a part of the initial 
Compton-produced photons to get into a cold matter and 
reflected. 

We thus fltted the XIS2, HXD-PIN, and HXD- 
GSO spectra simultaneously by a wabs* (diskbb+compPS) 
model, in the same manner as before. At this stage, 
the reflection option in compPS was suppressed, and the 
Gaussian component was not included. Here and here- 
after, we assume that the Compton cloud has a spherical 
geometry (parametcr=4) with an electron temperature To 
and an optical depth r, and that the seed photons to 
the compPS component are supplied by a cool standard 
disk having the same Tin (but a separate normalization) 
as the diskbb component. The implied picture is that 
some fraction of photons from a standard disk is fed to 
the Comptonizing corona, while the rest is directly visible 
as the spectral soft excess. This picture is supported by 
the presence of a tight correlation between the continuum 
slope and the reflection strength in BHBs (Zdziarski et al. 
1999; Gilfanov et al. 1999; Revnivtsev et al. 2001) and 
in AGNs, which indicates that the optically-thick disk, 
producing the reflection, is also causally related to the 
Comptonized emission. 




J , ■ , ■ < 
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Fig. 4. Simultaneous fittings to the time-averaged XIS2, 
HXD-PN, and HXD-GSO spectra of Cyg X-1, employing ther- 
mal Comptonization models. The derived best-fit parameters 
are given in table 1. (a) The same spectra as in figure 3a, com- 
pared with the best-fit wabs*(diskbb+conipPS+compPS+gau) 
model. The two compPS components are constrained to have 
the same Te, and a common set of refiection parameters, (b) 
Fit residuals from a model incorporating a single compPS con- 
tinuum without refiection, together with a diskbb. (c) The 
same as panel b, but refiection is incorporated, and a Gaussian 
for the Fe-K line is added, (d) The residuals correspond to 
the fit shown in panel a. Compared to panel c, the second 
compPS component is added. 

As shown in figure 4b, this model gave a poor fit with 
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'X^ jv = 3.9 for V = 357. Since figure 4b reveals the reflec- 
tion hump at ^ 30 keV, we activated the reflection option 
in compPS. This has nearly the same effect as the pexrav 
model, but with two additional free parameters; the ion- 
ization parameter ^ in units of erg cm s~^, and the inner- 
most reflector radius i?™^. The latter controls relativistic 
blurring of spectral features, and is expressed in units of 
the gravitational radius i?g = GAfeH/c^, with Mbh the 
BH mass and G the gravitational constant. The outer ra- 
dius was fixed at 10^i?g and the radial emissivity profile 
was assumed to depend on the radius r as (Fabian 
et al. 1989). The inclination angle, abundance, and tem- 
perature of the reflector were also fixed at 45°, 1 solar, 
and 10^ K, respectively. We retained the Gaussian com- 
ponent, because the compPS reflection does not include 
any emission line. Then, as shown in figure 4c, the fit was 
much improved {x^ /v = 2.6 with v ~ 351); the refiection 
hump is significantly present with ^ 0.2 and ^ ~ 0, 

although R\'^ was not constrained in this and all subse- 
quent fits. The fit goodness is comparable to that ob- 
tained with the wabs* (diskbb+cutof f pl+pexrav+gau) 
model (figure 3d). While the conventional model failed to 
constrain T^ut, the present model gave an electron tem- 
perature of To = 115 keV. 

In figure 4c, the fit still leaves a prominent data deficit 
in the highest XIS range. In addition, the XIS vs. HXD 
normalization ratio was obtained at 0.108, which exceeds 
the nominal value of ^ 0.092 (subsection 2.2); this is too 
large, even considering uncertainties in the correction fac- 
tors involved in the XIS and HXD data processing (sec- 
tion 2). We hence suspect that the actual source contin- 
uum is rising from ~ 10 kcV toward lower energies, with 
a steeper slope than indicated by the model, so that the 
fitting algorithm tried to absorb this mismatch by arti- 
ficially raising the model prediction to the XIS data. In 
other words, the broad-band continuum has a more con- 
cave shape over 2-20 keV (see figure 3b), than would be 
explained by the compPS component with reflection. This 
effect was previously noticed by Di Salvo et al. (2001) as 
a soft excess that cannot be accounted for by the cool 
disk emission. The same effect was noticed in the Suzaku 
data of GRO J1655-40 (Paper I), and led Takahashi 
et al. (2007a) to incorporate another compPS component. 
Such a "double-Compton" modeling was previously ap- 
plied successfully to the LHS of Cyg X-1 in several works; 
by Gierlinski et al. (1997) and Ibragimov et al. (2005) us- 
ing simultaneous Ginga, RXTE, and CGOR/OSSE data, 
and by Frontera et al. (2001a) to the BcppoSAX data. 

We hence tried a model incorporating two Comptonized 
continua, which are constrained to have the same seed- 
photon temperature (which in turn is the same as that 
of diskbb), the same hot-electron temperature, and the 
same reflection parameters, but are allowed to differ in 
their normalization and optical depth. The model is hence 
wabs* (diskbb+compPS+compPS+gau) . As shown in fig- 
ure 4d, the fit has been improved drastically to /v— 1.13 
(jy = 349); the XIS, PIN, and GSO spectra are reproduced 
with /v = 1.3, 0.9, and 0.9, respectively. The model 
has removed the data deficit in the highest XIS range, 



and reduced a residual structure near 2 kcV where the 
two compPS components now cross over. Furthermore, the 
XIS vs. HXD normalization has became 0.088, which is 
close to the value of 0.092 calibrated using SS Cyg (sub- 
section 2.2). We are therefore confident that the Suzaku 
data require two (or possibly more) Compton components. 
Although the fit is not yet formally acceptable, we regard 
it as satisfactory, because the data vs. model discrepancy, 
typically within 4% over the entire 0.7-400 keV range, 
is comparable to those obtained when we fit the Suzaku 
spectra of bright objects, such as the Crab spectra, with 
simple models. 

Our final model obtained in this way is shown in fig- 
ure 4a in the convolved form, and in figure 5a in the inci- 
dent vFv form. Its model parameters are given in table 1. 
Hereafter, we call the compPS components with the larger 
and smaller j/-parameters compPSh and compPSs, respec- 
tively, with the suffix "h" standing for "hard" and "s" for 
"soft". In figure 5, the compPSs component includes not 
only the scattered photons [ex 1 — exp(— r)], but also those 
seed photons which traversed the Compton cloud without 
scattered [(xexp(— r)]. 



■I I I 




tUb) GRO J 165 5^40 




Energy (keV) 

Fig. 5. (a) The inferred best-fit vFv spectrum of Cyg X-1, 
shown with the absorption removed. It corresponds to panels 
a and c of figure 4, with the model parameters detailed in 
the column labeled "Average" of table 1. Red, blue, orange, 
green, and purple specify compPSh, compPSs, diskbb, reflection 
components, and a Gaussian for the Fe-K line, respectively, 
(b) The same as panel a, but for GRO J1655— 40 observed 
with Suzaku on 2005 September 22 and 23. It is identical to 
figure 7 of Paper I, except for the removal of the absorption. 
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3.3. Examination of the Comptonization model results 

Although the analysis conducted in the preceding sub- 
section has provided a plausible representation of the 
time-averaged spectra, different spectral models could of- 
ten degenerate. Therefore, let us briefly examine our final 
spectral model for its implications and physical consis- 
tency, leaving detailed discussion to section 5. 

The obtained absorption, (6.6^0 3) x 10^^ cm^^, is close 
to the value of (5.3 ± 0.2) x 10^^ cm~^ derived by Dotani 
et al. (1997) with ASCA. The compPSh and compPSs com- 
ponents, crossing at ~ 4 keV, are characterized by photon 
indices of ^1.6 and ^ 2.4, respectively, with the former 
consistent with the approximate X-ray slope of Cyg X- 
1 measured in energies of a few to a few tens keV. The 
underlying parameters, To ~ 100 kcV, r of order unity, 
and Q/2n ^ 0.4, are consistent with a large number of 
past measurements (e.g., Gierlinski et al. 1997; Di Salvo 
et al. 2001; Frontera et al. 2001a; Zdziarski & Gierlinski 
2004; Ibragimov et al. 2005). 

The diskbb component has a raw innermost radius as 
Till = 180^50° km, which yields the best-estimate actual 
inner radius of i?in = 250^gj!)'' km; this is obtained by ap- 
plying the nominal correction factors in Makishima et al. 
(2000) (including a color hardening factor of 1.7; Shimura 
& Takahara 1995) to r-m, and correcting for the inclina- 
tion assuming i ~ 45° , as 

Rin = l.lSrin/vcosz . (1) 

This i?in is larger than the value of i?in ^ 100 km mea- 
sured in the high/soft state with ASCA (Dotani ct al. 
1997; Makishima et al. 2000), suggesting a disk trunca- 
tion at a larger radius than in the high/soft state (see 
also subsection 5.4). The innermost disk temperature, 
0.19 keV, is a factor 2.3 lower than that measured in the 
high/soft state (Dotani et al. 1997), consistent with a re- 
duced accretion rate and possibly a larger innermost disk 
radius. These consistencies support (though not necessar- 
ily prove) our basic assumption that the seed photons for 
Comptonization are provided by the optically-thick disk, 
as proposed previously (Zdziarski et al. 1999; Gilfanov 
et al. 1999; Revnivtsev et al. 2001). Although the bolo- 
metric luminosity of the directly-visible disk emission is 
I9tl'^% of the overall 0.7-300 keV luminosity, the energy 
generation in the cool disk is approximately twice higher, 
when the seed photon luminosity is added. Then, the en- 
ergy generation in the the cool disk becomes comparable 
to that in the hot corona. 

While the compPSh component carries the dominant 
hard continuum, compPSs successfully explains the con- 
tinuum steepening in energies below a few keV, or equiv- 
alently. the soft excess that is too hot to be explained 
by the cool disk emission (Di Salvo et al. 2001). These 
two Compton continua are inferred to have optical depths 
which differ by a factor of 3.9 ±0.5. Their normalizations, 
originally expressed as the disk area A supplying the seed 
photons, have been converted in table 1 to its innermost 
radius Rf°°'^ as 



^ = 1.18\/^ . (2) 

Though similar to Rm of equation (1), i??™'^ is not 
corrected for the inclination, since we presume the 
Comptonized emission to be approximately isotropic. As 
seen in table 1, the value of i?f°°'^ associated with compPSs 
is comparable to i?in, while that of compPSh is consider- 
ably smaller. We hence obtain an implication that ~ 7 
times more seed photons are fed to compPSs than to 
compPSh: a further discussion continues in subsection 5.4. 

Although the final model is generally successful, the 
obtained XIS vs. HXD normalization, 0.088, is ^ 5% too 
small. When this ratio is fixed to the calibration value 
of 0.092, the PIN data exhibit some positive and negative 
residuals at ~ 15 keV and ^ 40 keV, respectively. These 
can be eliminated by allowing compPSs and compPSh to 
have a larger (~ 1.0) and a slightly smaller (^ 0.3) value 
of ri/27r, respectively, with the raw chi-square reducing by 
13 (for ly = 349). This is because the reflection associated 
with compPSs appears at lower energies than those with 
compPSh (figure 5a). There is thus some indication that 
the softer Compton signals are more strongly refiected by 
cold materials than the harder one. 

When the two compPS components are allowed to take 
different values of Tc, tolerances of the model parameters 
increase, but the fit goodness remains the same. This is 
because the high-energy roll over of compPSs is masked 
by compPSh. Therefore, we do not adopt this modeling, 
although the data do not necessarily exclude the presence 
of multiple electron temperatures. 

4. Analysis of Spectral Changes Associated with 
the Random Variation 

4.1. Methods 

Now that the time-averaged spectra have been quanti- 
fied in section 3, we proceed to the latter half of our data 
analysis, namely attempts to understand how the spec- 
trum changes as the source varies randomly; a similar at- 
tempt was carried out previously by Takahashi, Inoue & 
Dotani (2001). For this purpose, it is essential to sort 
the broad-band spectra, particularly those of the HXD, 
into two subsets, according to whether the instantaneous 
source count rate is higher or lower than a certain mean 
value. However, the HXD has a relatively small effec- 
tive area 230 cm^ at 100 keV; Takahashi et al. 2007c), 
so the average count rates of Cyg X-1 are ~ 47 c s~^ 
and - 20 c s-Mn the 10-70 keV (with PIN) and 70-400 
keV (with GSO) bands, respectively (subsection 2.1; fig- 
ure lc,d). Then, even in the PIN range where the NXB is 
relatively minor, Poisson fluctuations of signal counts with 
1 s integration would become comparable to the intrin- 
sic variations of Cyg X-1 (~ ±20% in root-mean-square 
amplitude). Making the integration time longer would 
not help, because Cyg X-1 varies with an approximately 
"white" power spectrum on time scales longer than a few 
seconds (e.g., Belloni & Hasinger 1990), and hence ampli- 
tudes of the intrinsic and Poisson variations will scale in 
a similar manner when the integration time is increased. 
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Fortunately, we can utilize the XIS for our purpose, 
because one XIS sensor received ^ 60 c s~^ even after the 
image center has been removed to avoid photon pile up 
(subsection 2.2; figure la,b), and the variation of Cyg X-1 
is known to be crudely synchronized across the soft X-ray 
to hard X-ray energies (Negoro et al. 1994). Since the 
present XIS data have a 1 s time resolution (subsection 
2.1), we can judge every second, referring to the XIS data, 
whether the source is brighter or fainter than a mean, and 
sort the XIS and HXD data into the required two subsets. 
In practice, we have constructed the following procedure 
to do this. 

1. In the same way as in section 2.4, a 0.3-10 keV 
XIS0+XIS3 light curve {CJ (with i the bin num- 
ber) is produced throughout the observation, with a 
bin width of tb > 1 s. 

2. The calculation of {Ci} excludes those 1-s time bins 
in which either XISO or XIS3 (or both) suffered 
telemetry saturation. 

3. Every T s (T > t,,), a mean XIS0+XIS3 count rate 
C is calculated from {Ci}. 

4. Two sets of good time intervals (GTIs), called "high- 
phase" (HP) and "low-phase" (LP), are defined, by 
judging every th s whether d is higher than C + na, 
or lower than C — na, respectively. Here, cr is 1- 
sigma Poisson error and n > is a discriminator 
level. 

5. Using the same two GTIs, the XIS2 data are ac- 
cumulated into two spectra, called HP (high phase) 
and LP (low phase) spectra. The data are discarded 
\i C - na < Ci < C + na . 

6. Employing the same two GTIs, the HXD-PIN data 
are sorted into HP and LP spectra; so are the HXD- 
GSO data. 

7. Using the two sets of GTIs, two background data 
sets are accumulated for HP and LP, and are 
subtracted from the corresponding on-source HXD 
spectra. 

Thus, the procedure involves three adjustable parame- 
ters, th, T, and n. Among them, T should be long enough 
to suppress Poisson errors, but short enough to compen- 
sate for background changes. The threshold n should also 
be optimized, because a larger n leads to a decrease of 
live time, while a smaller n means that the data sorting 
is diluted by the Poisson fluctuation in {Ci} which is not 
related to the intrinsic source variation. Here, for simplic- 
ity, we have chosen = 1 s (the highest time resolution 
of the XIS); T = 200 s, and n = as our baseline condi- 
tions. This is because a larger value of n, such as n = 1, 
was found to enhances the HP vs. LP ratio, but also to 
degrade the data statistics, thus hampering us to more ac- 
curately separate the HP and LP model parameters. Since 
we have chosen to use the XIS0-I-XIS3 data for the inten- 
sity judgement, we derive the intensity-sorted XIS spectra 
only from the other FI sensor, namely XIS2. This divi- 
sion of task is needed because the intensity-sorted spectra 
from XISO and XIS3 would include those time bins when 
their summed counts are higher (or lower) due simply to 



statistical fluctuations, which amount to ^ 8% (1 sigma) 
at the 1-s integration. The intensity-sorted spectra from 
XISO and XIS3, if to be utilized in the spectral fitting, 
would need to be corrected for this effect. The XIS2 data, 
in contrast, are free from this problem, because Poisson 
errors from different detectors are independent. 

4-2. Auto- correlation and cross- correlation functions 

Before actually conducting the analysis described in 
section 4.1, we briefly perform auto-correlation function 
(ACF) and cross-correlation function (CCF) analyses, in 
order to examine whether the GTIs sorted in reference to 
the energy band below 10 keV is really applicable to the 
extraction of varying spectral components in harder ener- 
gies. We utilize raw light curves with neither background 
subtraction nor dead-time correction. 

Panels (a) and (b) of flgure 6 show ACFs and CCFs 
from the present Suzaku data, respectively, calculated em- 
ploying a bin width of 1 s. They were calculated using 
data intervals each 256 s long, and then averaged over 
ensemble. The XIS telemetry saturation was treated in 
the same way as item 2 of the previous subsection (using 
tb = 1 s). Thus, the XIS vs. PIN and XIS vs. GSO CGFs 
both exhibit clear peaks at zero time lag, with correla- 
tion coefficients of ~ 0.75. This ensures that the intrinsic 
variation on a time scale of 1 s is well correlated over the 
entire energy range in which we are conducting our anal- 
ysis, and validate our method described in the previous 
subsection. 

Although the XIS data do not have time resolution 
higher than 1 s, the HXD data afford it down to 61 fis 
(Takahashi et al. 2007c; Terada et al. 2007). We calcu- 
lated ACFs and CCFs from the PIN (10-60 kcV) and 
GSO (60-200 keV) light curves accordingly, with 0.1 s 
resolution and 409.6 s data intervals. The results, shown 
in panels (c) and (d) of figure 6, reveal that the varia- 
tions are detected on shorter time scales, where the PIN 
vs GSO correlation is still peaked at zero-lag. The varia- 
tions in the two energy bands are again well synchronized, 
with the correlation coefficient of 1. 

In addition to the basic properties described above, 
these ACFs and CCFs reveal one interesting property: 
namely, the left-to-right asymmetry in the CCFs, seen 
on both time scales, in such a way that the correlation 
is stronger toward hard-lag sense. This reconfirms the 
previous report of the same phenomenon in early days 
by Nolan et al. (1981), then using Ginga by Miyamoto 
et al. (1991), and further with RXTE (e.g., Nowak et al. 
1999a; Maccarone et al. 2000). On the other hand, the 
present results do not clearly reconform another property, 
reported previously (e.g., Negoro et al. 1994; Feng, Li & 
Chen 1999; Maccarone et al. 2000), that ACFs in softer 
energies have wider correlation peaks, and hence longer 
correlation lengths. We leave this issue to future Suzaku 
observations using higher XIS time resolution (employing 
so-called PSUM mode); then, we will also be able to study 
how the spectrum differs between the rising and falling 
phases of "shots", as suggested previously (Negoro et al. 
1994). 




Fig. 6. Auto-corrclation functions (ACFs) and cross-correlation functions (CCFs) of the present Suzaku data. The effects of Poisson 
errors were removed, and the CCFs were normahzed so that the value at zero lag directly represent the correlation coefficient. Positive 
time differences in CCFs mean that softer signals lead the harder ones, (a) ACFs from the X1S2 (0.7-10 keV; black), HXD-PIN 
(10-60 keV; red), and HXD-GSO (60-200 keV; green) data, calculated using 1-s binned light curves, (b) CCFs of the PIN (red) 
and GSO (green) data, calculated against the 0.7-10 keV XIS2 light curve in the same way as panel a. (c) ACFs calculated using 
0.1-s binned PIN (10-60 keV; red) and GSO (60-200 keV; green) light curves, (d) CCFs of the GSO light curve against that of PIN, 
calculated in the same way as panel c. 




Time (sec) 

Fig. 7. Background- inclusive XIS (top: 0.7-10 keV), HXD-PIN (middle: 10-60 keV), and HXD-GSO (bottm: 60-200 keV) light 
curves of Cyg X-1 with 1-s binning, sorted according to the instantaneous XIS0-I-XIS3 counts. Red and blue data bins represent 
those when the XIS0-I-XIS3 counts, d, arc higher and lower than the local 200 s average, C, respectively. No correction is made for 
the instrumental dead times. 
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^.3. Intensity- sorted spectra 

Now that the examination of CCFs have vahdated the 
method presented in section 4.1, we have appHed it to the 
present Suzaku data. Figure 7 gives an example of 1-s 
bin hght curves of Cyg X-1, sorted according to the above 
method using the basehne condition; ib = 1 s, T = 200 s, 
and n = (section 4.1). In this particular example, the 
average XIS count rate used as the HP vs. LP threshold is 
140 and 148 ct s~^ in the former and latter 200 s, respec- 
tively. These are ^ 15% higher than the XIS2 counting 
rate indicated by figure 1 (panels a, b), because the XIS2 
optical axis is more offset from the 1/8-window centroid 
than those of XISO and XIS3, causing a larger photon loss. 
The XIS light curve reveals the well known flickering be- 
havior of Cyg X-1, while those of the HXD suffer larger 
Poisson errors so that the intensity sorting on a time scale 
of 1 s is impossible without the XIS data. 

Figure 8 shows the HP and LP spectra derived in this 
way. There, we subtracted the corresponding backgrounds 
accumulated using the respective GTIs (item 7 in sec- 
tion 4.1): the utilized two background sets agree within 
their statistical errors. While the spectra are presented 
without dead time correction, we correct the model noma- 
lizations (in table 1 and elsewhere in text) for the time- 
averaged dead time of 7%. Although the HP and LP 
spectra should have different dead times, the difference 
is calculated to be within ~ 2%, because the dead time is 
mainly caused by background counts. This difference is 
expected to reduce the HP spectral normalization by at 
most 2%, keeping the spectral shape intact. 

Figure 8b gives the ratios between the HP and LP spec- 
tra, with several important implications. First, the aver- 
age ratio becomes ^ 1.3, in agreement with the typical 
variation amplitude of Cyg X-1. Second, the ratio in- 
creases at the highest XIS energies band, but this is due 
to still remaining photon pile-up effects in the XIS when 
the source gets brighter, and should be ignored. Third, 
the spectrum gets slightly softer when the source bright- 
ens up, as evidenced by a gradual decrease of the HP/LP 
ratio toward higher energies; this qualitatively agrees with 
previous results from, e.g., Ginga (Negoro et al. 1994). 
Finally, the HP /LP ratio exhibits a rather abrupt decrease 
at 1 keV toward lower energies. 

In order to see the effects of changing the selection con- 
dition, we repeated the intensity sorting using ib = 3 s 
and 8 s. The results obtained with tb = 3 s were essen- 
tially the same as those with the baseline case {t\, = 1 s). 
When employing ib = 8 s, the HP/LP ratio was reduced 
to ^ 1.2 on average, although its energy dependence was 
very similar to the baseline case. From these, we retain 
ib = 1 s. We also examined different Poisson thresholds, 
n. When rt = (the baselline case), 1, and 2 are used, the 
HP GTIs yielded an overall exposure of 7.3, 5.2, and 3.3 
ks, respectively (with LP GTI being comparable), while 
the energy-averaged HP/LP ratio turned out to be ^ 1.3, 
~ 1.4, and ~ 1.5, respectively. We hence retain ri = 0, since 
it maximizes the signal-to-noise ratio of the HP minus LP 
difference spectra. 
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Fig. 8. (a) The HP (high phase; black) and LP (low phase; 
red) spectra of Cyg X-1 recorded by XIS2, HXD-PIN, and 
HXD-GSO, while XIS0-I-XIS3 data are used in the intensity 
judgement. The corresponding HXD backgrounds have been 
subtracted, (b) The ratios between the HP and LP spectra. 

Returning to the HP and LP spectra obtained with the 
baseline condition, we tried grossly quantifying them using 
the conventional model employed in section 3.1. To avoid 
the complexity in softer energies which arisses from the 
iron line, the soft excess, and the softer Compnotization 
continuum, we fitted only the HXD (PIN and GSO) spec- 
tra, using the cutoff pl+pexrav model. The model re- 
produced the HP and LP spectra successfully both with 
x/v ^ 1, and yielded the consistent reflection strength; 
rj/27r = 0.16t!^:g^ in HP and 0.16 ±0.06 in LP. (These are 
smaller than the value of ~0.4 obtained with the time- 
averaged spectra using the Compnonization models, be- 
cause a cutoff power-law model is more convex than an 
unsaturated Compton continuum.) 

To highlight the difference between the two spectra, we 
next fitted them jointly, with VI/2tt fixed at 0.16. As show 
in table 2 (column under an entry "Free"), we were un- 
able to find statistically significant differences either in F 
or E'cut- Therefore, we further constrained the HP and LP 
spectra to have either the same F or the same i?cut- Then, 
as shown in the last two columns of table 2, the join fit 
remained acceptable without preference between the two 
cases, and the unconstrained parameter showed a statis- 
tically significant difference between the two spectra. We 
therefore conclude that the HP and LP spectral have dif- 
ferent shapes as suggested by their ratio (figure 8b), and 
that either i?cut or F (or both) is changing, although we 
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cannot determine which is more variable. 
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Fig. 9. Background-subtracted and dead-time uncorrected 
HP (panel a) and LP (panel b) spectra of Cyg X-1, fitted with 
the wabs* (diskbb+compPS+compPS+gau) model in the same 
way as figure 4. Panel c is for the difference spectra between 
HP and LP. 



are found to generally bracket that from the average spec- 
trum. The inferred best-fit models, in the i^Fiy form, are 
given in figure 10. 
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4-.4- Changes of the Compton parameters 

Following the approximate quantification of the 
intensity-sorted HXD spectra (subsection 4.3), we fitted 
the full (XIS2, PIN, and GSO) HP and LP spectra with 
the wabs*(diskbb +compPS+compPS+gau) model which 
has been found to be the best representation of the time- 
averaged spectrum (section 3.2). In the fitting, we fixed 
A^H = 6.6 X 10^^ cm~^ as obtained with the average spec- 
tra. In reference to the preliminary fittings conducted in 
the previous subsection, we also fixed all the reflection pa- 
rameters to the values obtained with the average spectra. 
Then, as presented in figure 9, the model has success- 
fully reproduced both spectra. The derived parameters 
are listed in table 1, in comparison with those from the 
averaged spectrum. Thus, the LP and HP values of each 
free parameter (except a few rather unconstrained ones) 



Fig. 10. The same as figure 5a, but obtained from the HP 
(panel a) and LP (panel b) spectra. Panel c is the fit to the 
HP minus LP diff'eronce spectra. 



Table 1 and figure 10 allow us to derive several im- 
portant inferences. First, let us consider the observed 
spectral softening from LP to HP (figure 8b). Although 
this effect might suggest an increased relative contribu- 
tion of compPSs in HP, this is not the case. In fact, as 
the source varies, the seed-disk radius of compPSs is kept 
at 3± 1 times that of compPSh (table 1), and their cross- 
over energy stays at ~ 4 keV (figure 10). In search for an 
alternative explanation, figure 11 compares the compPSh 
parameters between HP and LP on the plane of Tc against 
the optical depth Th. Thus, the HP and LP spectra be- 
come distinct on this plane, with implication similar to 
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what was derived in subsection 4.3; the two phases can- 
not be represented by a common set of Tc and Th. If, 
e.g., To is the same, then Th must be lower in HP than in 
LP; if instead Th is assumed to be the same, To must be 
higher in LP. The two phases are more clearly (though not 
completely) separated in j/, because it can be accurately 
specified by the observed spectral slopes. We may, for sim- 
plicity, express that the hard continuum is "less strongly 
Comptonized" in HP than in LP. This is an important re- 
sult enabled by the broad-band capability of Suzaku. As 
to compPSs, we cannot find significant differences between 
LP and HP beyond their relatively large errors. 
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Fig. 11. Confidence contours of the compPSh paramters in 
HP (solid lines) and LP (dashed ones), shown on the plane 
of Te versus the optical depth rjj. Red and green contours 
represent 90% and 99% confidence levels, respectively. Blue 
lines indicate constant values of y. 

Next, the source variation has been found to approxi- 
mately conserve the innermost disk temperature at Tin = 
0.19 keV (table 1). In the present modeling, this quantity 
is assumed to be common to the directly visible diskbb 
component, and the seed photon sources for the two 
compPS components. Therefore, it is determined mainly 
by local shapes of the spectral soft excess (figure 10). 
However, the broad-band capability is again very impor- 
tant, because it allows us to accurately fix the Compton 
continua and the reficction hump, and hence to unam- 
biguously model the soft excess. 

As a third result from this analysis, we can at- 
tribute the LP to HP intensity increase to an increase 
in the spherically-averaged seed-photon luminosity, which 
is given as 

isocd = 47rcro'7lti|i?in'l1i + Rtn^-s } i (3) 

here, ctq is the Stefan-Boltzmann constant, while i?f,j'i'|, 
and Rf^^'^g respectively represent the seed-disk radius of 
compPSh and compPSs, as defined by equation (2). When 
Tseed is estimated table 1 via equation 3, its nominal error 
ranges overlap between HP and LP, due to large errors in 
Rf°°'^. However, in reality, i?;^'^'^ is tightly anti-corrclated 
with Tin so as to preserve the seed-photon fiux (for compPS 
in particular). Therefore, the combined error range of 



Tsood becomes much smaller. A detailed confidence-range 
examination yields Tgcod — (S.S^Qg) x lO'^^ erg s~^ in HP, 
and {6.9t°-l) x lO^^ erg s"^ in LP. 

Finally, the bolometric luminosity of the directly- visible 
diskbb component, calculated after equation (3) as 

iraw = 4^aoIlti?in' (4) 

and presented in table 1, is consistent with being the same 
between LP and HP, within rather large errors which 
arise mainly via a strong coupling between diskbb and 
compPSs. These error ranges associated with Traw, which 
also take into account the temperature vs. radius anti cor- 
relation, are such that a significant decrease (by as much 
as a factor of ~ 5) in Traw from LP to HP is allowed, while 
its increases should not exceed ~ 20%. We therefore infer, 
as our 4th result, that the directly visible disk luminos- 
ity is likely to be less variable than the hard X-ray flux, 
and there remains a possibility that it even anti-correlates 
with the hard X-ray flux. 

4-5. Difference spectra 

To reconfirm and reinforce the results obtained in 
the preceding subsection, we directly subtracted the raw 
(i.e., background- inclusive and dead-time uncorrected) LP 
spectra from those in HP. Since the source resides in ei- 
ther phase for a typical duration of several seconds which 
is much shorter than the time scale of the NXB variation, 
we expect the NXB to accurately cancel out by this sub- 
traction, making "difference spectra" free from the NXB 
modeling errors. In fact, over the highest energy of 530- 
600 keV where the signal is considered negligible (< 1% of 
the NXB), the HP minus LP count rate of HXD-GSO has 
been found to be (— 0.6± 1.0)% of the time-averaged mean 
rate. This also confirms that the potentially higher (by 
at most 2%; subsection 4.3) dead time in HP than in LP 
can be neglected. The derived difference spectra arc pre- 
sented in figure 9c. As expected, they have several times 
smaller normalizations than the HP and LP spectra, and 
are generally softer. 

As before, we fitted the difference spectra using the 
progressively complex Comptonization models, but fixing 
A^H = 6.6 X 10^2 cm-2, Tin = 0.19 keV, fi/27r = 0.4, ^ = 0, 
and ~ 2500i?g. The systematic errors were not in- 
corporated due to lower statistics. Then, a single compPS 
continuum, together with reflection and the diskbb com- 
ponent, gave a poor fit with jv^ 3.36 for v = 67. The fit 
goodness was improve to ^ jv = 3.18 iy = 66) by activat- 
ing the Gaussian component, of which the center energy 
and the width are fixed at 6.3 keV and cr = 1.0 kcV, respec- 
tively. Finally, by incorporating the second compPS com- 
ponent, the fit has become fully acceptable {x^ /v = 1.05 
for V = 64), yielding the parameters given in table 1 (col- 
umn before the last). Thus, the model has again the same 
construction as those for the time-averaged, HP, and LP 
spectra. The best-fit model is presented in figure 9c in 
the convolved form, and in figure 10c in the incident vFv 
form. Below, we examine the fit results. 

In the difference spectra, the compPSs component has 
been found to have again ~ 3 times larger seed-disk radius 
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than compPSh. Since this ratio is close to those found in 
the average, HP, and LP spectra, we reconfirm that the 
soft-photon inputs to the two Compton clouds are vary- 
ing roughly in proportion to each other. In contrast, the 
value of To ~ 60 keV is considerably lower than was found 
before. One possible scenario from these results is that 
the Compton-cloud temperature in fact has a distribu- 
tion, which becomes more weighted toward lower values 
as the source flares up, so that the difference can be ap- 
proximated by a single lower temperature. 

In subsection 4.3, we confirmed (though using the 
conventional model) that the HP and LP spectra have 
roughly the same ri/27r. This is reinforced by the re- 
sults obtained subsequently, namely the successful double- 
Compton modelings of the HP, LP, and the difference 
spectra, all with r2/27r fixed at 0.4. Therefore, the reflec- 
tion component is thought to be fully catching up with the 
continuum variations on the 1-s time scale. Also, the iron- 
line photon flux is varying clearly on this time scale; its 
variation is consistent with being directly proportional to 
the continuum flux, because all the four spectra of Cyg X- 
1 presented in table 1 have, within errors, the same Fe-K 
line equivalent width of ^ 300 eV . These results are con- 
sistent with Rcvnivtscv ct al. (1999), who reported that 
the reflection features (Fc-K line and smeared edge) start 
failing to follow continuum variations above several Hz in 
frequency. 

Although the flt to the difference spectrum incorporates 
the diskbb component, the final model given in table 1 
does not necessarily require it: only an upper limit can 
be imposed on Rin, because the softest end of the differ- 
ence XIS spectrum can be accounted for by compPSs as 
well. We in fact confirmed that a negative normalization 
of diskbb (i.e., a decrease of this component from LP to 
HP) is allowed up to ^ 40% (in absolute value) of those 
in LP. This is consistent with the results obtained in the 
previous subsection. 

4-6. Summary of spectral changes 

Through the analyses in the preceding three subsec- 
tions, the spectra changes associated with the fast (1 to 
200 s) intensity variations have been quantified in the fol- 
lowing manner. 

1. The intensity increase from LP to HP is explained 
by an increase in the seed-photon luminosity, igccd- 

2. The spectral softening from LP to HP can be at- 
tributed to weakening in the Comptonization of 
compPSh (either small r, or lower Te, or both). 

3. Unlike Lseed, the directly visible disk luminosity, 
Lraw) does not increase significantly in HP, and 
could even decrease. 

4. As the source varies, the innermost temperature of 
the cool disk is kept approximately constant at Tin = 
0.19 keV. 

5. On these time scales, the Fe-K line and the reflec- 
tion hump follow the continuum variations, because 
il/27r and the Fe-K line EW are both consistent with 
being constant. 



5. Discussion 

5.1. Summary of the Data Analysis 

The Suzaku observation of Cyg X-1, conducted for 
an exposure of 17 ks at a 0.7-300 keV luminosity of 
4.6 X 10^^ erg s~^, has provided one of the highest-quality 
LHS spectra of this prototypical BHB. The XIS, HXD- 
PIN, and HXD-GSO spectra, altogether covering an ex- 
tremely broad energy band of 0.7-400 keV, have been ex- 
plained simultaneously and consistently, by invoking two 
basic constituents of the X-ray emission. One is the hot 
Comptonizing corona with Tc ~ 100 keV, which manifests 
itself in the high-energy spectral roll over detected with a 
high signiflcance by HXD-GSO. The other is the optically- 
thick cool disk with Tin ^ 0.2 keV, which produces the soft 
excess and the Fe-K line in the XIS band, as well as the 
reflection hump in the HXD-PIN band with f7/27r~ 0.4. 

The data are fully consistent with the cool disk sup- 
plying seed photons to the Compton cloud, although 
other possibilities may not be necessarily excluded. In 
the present modeling, the soft photons from the cool 
disk are considered to reach us either directly as the 
soft excess (the diskbb component), or through weaker 
Comptonization with y ~ 0.3 (compPSs), or after experi- 
encing stronger Comptonization with y ^ 1.15 (compPSh). 

Assembling together the two Compton components 
(compPSh, compPSs) and the three cool-disk related ones 
(diskbb, reflection, and gaussian), we have successfully 
reproduced the overall continuum shape, including the 
complicated Fe-K line and soft-excess regions. The Fe- 
K line is inferred to be broadened only weakly to ~ 1 keV 
(in Gaussian a), without any evidence for relativistically 
broadened wings of which the detection is claimed from 
some BHBs (Miller 2006). As seen in flgure 5, the in- 
ferred incident spectrum shows too complicated a shape, 
all over the energy range, to be approximated by a single 
power-law. 

Analyzing in section 4 the characteristic flare-up be- 
havior of Cyg X-1 on a time scale of ~ 1 s, we found 
that the wabs* (diskbb+compPS+compPS+gau) modeling 
applies also to the HP, LP, and the difference spectra. 
The flare-up behavior can be understood primarily as an 
increase in the seed-photon supply, accompanied by a de- 
creased Comptonization. While the seed-disk luminosity 
thus increases in HP, the luminosity of the directly visi- 
ble diskbb does not increase significantly. As the source 
varies. Tin of the cool disk is kept approximately constant. 

5.2. Comparison with the BeppoSAX results 

Applying the double-compPS modeling to the 
BeppoSAX data of Cyg X-1 in the LHS covering up to 
^ 200 keV (whereas Suzaku reached 400 keV), Frontcra 
et al. (2001a) derived Te = 59±5 keV, y = 0.89 ±0.01, 
and ri/27r = 0.25 ± 0.04 for compPSh (their Compton 
component 1), and To = 42 ± 19 keV and y = 0.15 ± 0.01 
for compPSs (their copmPS component 2). Except for 
some differences in the model parameters, the Suzaku 
and BeppoSAX results are hence consistent in that 
the broad-band spectra can be reproduced by a pair of 
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Comptonized components. The two results grossly agree 
on the soft-excess parameters as well, although Frontera 
et al. (2001a) used a blackbody to model the soft excess 
and the sccd-photon source. 

The BeppoSAX observation was conducted at a 0.5-200 
keV unabsorbed flux of 4.2 x 10~^ erg cm~^ s~^, or a lumi- 
nosity of 3.1 X 10'^'^ erg s^^ (in the same band at a distance 
of 3.2 kpc). The corresponding values from the present 
Suzaku observation are higher, 5.9 x 10~^ erg cm~^ s~^ 
(4.4 X 10'^^ erg s^^). Then, we would expect Suzaku to 
observe a lower To, based on the general understanding 
that Te decreases toward higher luminosities due to en- 
hanced Compton cooling (e.g., Esin ct al. 1998; Yamaoka 
et al. 2005; Zdziarski & Gierhnski 2004). Nevertheless, 
we in fact measured a 1.6 times higher Tq than BeppoSAX. 
One possible explanation to this behavior is that Cyg X-1 
may not have fully returned to the LHS at the time of the 
BeppoSAX observation, which was conducted soon after 
the source transition from the high/soft state. 

Prior to the present work, the double-Compton model 
was applied successfully to wide-band LHS spectra of 
Cyg X-1 (Gierhnski et al. 1997; Frontera et al. 2001a; 
Ibragimov et al. 2005; though some explicitly invoking 
multiple values of To) and GRO J1655— 40 (Paper I, sub- 
section 5.3). Therefore, we regard this model as a promis- 
ing description of the LHS spectra of BHBs. As discussed 
in Paper I, the mixture of two Compton optical depths 
(plus the direct diskbb component) allows either a "spa- 
tial" or a "time-domain" interpretation. In the former 
case, some fraction of the cool disk is considered directly 
visible, while the remaining part is covered by a hot elec- 
tron cloud that has thinner and thicker portions. The lat- 
ter (time domain) alternative assumes that the Compton 
cloud is rapidly varying among three typical conditions, 
represented by r = 0, t 0.4, and r ~ 1.5. The ac- 
tual condition may be a combination of these two sim- 
plified cases. Of course, as pointed out previously (e.g.. 
Paper I), the double-Compton modeling may be an ap- 
proximation to more complex conditions involving more 
than two Compton optical depths, or even a continuous 
distribution in r as suggested by some models (e.g., Liu 
et al. 2002). With the present data, it is virtually im- 
possible to distnguish which is more likely, a continuous 
distributin in t or just two optical depths. 

5.3. Comparison with GRO ,11655—^0 

A comparison between Cyg X-1 (with i ~ 45°) and 
GRO J1655-40 (with i - 70°) would help elucidating 
the source geometry (section 1). After a brief attempt 
made in Paper 1 using the BeppoSAX results on Cyg X- 
1 (Frontera et al. 2001a), here we perform a more accu- 
rate comparison between the two objects using the Suzaku 
data on both. For this purpose, figure 5b reproduces 
the best-fit wabs* (diskbb+compPS+compPS+gau) model 
for GRO J1655— 40, deduced in Paper I using the Suzaku 
data. On that occasion (2005 September 22 and 23), 
GRO J1655-40 was in the LHS with a 0.7-300 keV lumi- 
nosity of 4.9 X 10^^ erg s~^ (Paper I), which is an order of 
magnutude lower than that of Cyg X-1 in the present ob- 



servation (both assuming isotropic emission). However, 
the difference becomes a factor of 3-5 when the lumi- 
nosity is normalized to the Eddington value, because the 
BH in Cyg X-1 has a mass of 12 — 20 Mq, while that in 
GRO J1655-40 has - 6.5 Mq (Paper I). 
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Fig. 12. Background-subtracted XIS (black), HXD-PIN 
(red), and HXD-GSO (green) spectra of GRO J1655-40 
(Paper I), divided by those of Cyg X-1 obtained in the present 
study. 

Although the two models in figure 5 generally look alike, 
they differ in details. To grasp such differences in a model- 
independent manner, we divided the XI2, HXD-PIN, and 
HXD-GSO spectra of GRO J1655-40 by the correspond- 
ing spectra of Cyg X-1, and obtained the results shown 
in figure 12. Below, we attempt to identify basic features 
in this spectral ratio, and to interpret them in terms of 
the spectral decomposition in figure 5 and table 1. Like in 
figure 8b, the behavior of the ratio at the highest XIS en- 
ergy range should be ignored, where the GRO J1655— 40 
data may suffer from photon pile- up effects (Paper I). 

The spectral ratio in figure 12 is roughly constant over 
'--^ 3 to ~ 30 keV, but decreases significantly below ^ 2 
keV. This effect is too large to be explained by their dif- 
ference in absorption, A^h = 6.6 x 10^^ cm~^ in Cyg X-1 
and 7.4 x lO^i cm^^ in GRO J1655-40. Instead, it can be 
attributed to the stronger soft excess of Cyg X-1 (figure 5), 
as the diskbb component of Cyg X-1 is an order of mag- 
nitude more prominent than that of GRO J1655— 40 when 
normalized to the overall continuum (while Tin is nearly 
the same). This in turn may be ascribed primarily to 
their difference in the cosi factor, cos(45°)/cos(70°) =2.1, 
because we expect the direct disk emission to scale as 
cx cosi, while the Comptonized emission to be approxi- 
mately isotropic. Strictly speaking, this estimate needs 
some caution. In subsection 4.4, we introduced two quan- 
tities, Traw and Tsood- Because i?in used to calculate 
Traw has already been corrected for the cosz factor while 
Rl'^'^ is not (subsection 3.2, table 1), we expect the ratio 
Liaw/Lsccd to be approximately inclination independent. 
Nevertheless, the ratio was measured to be ^ 0.44 in Cyg 
X-1, while < 0.19 in GRO J1655-40 (table 1). That is, the 
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soft excess in GRO J1655— 40 is even weaker than would 
be predicted by the assumed i = 70°. This is discussed 
later in subsection 5.5. 

In figure 12, we observe a small dip at ~ 6.5 kcV. 
Evidently, this feature arises because the Cyg X-1 spec- 
trum exhibits a stronger Fe-K line (with an equivalent 
width of ~ 300 eV) than that of GRO J1655-40 (- 80 
eV). This property is again ascribed, at least qualitatively, 
to the inclination effects: the cool flat disk will intercept 
a fraction of the Comptonized photons, and produce the 
fluorescent Fc-K photons, of which the visibility must be 
proportional to cosz. Quantitatively, however, the dif- 
ference in the equivalent width between the two sources 
amounts to a factor of ^ 4, which exceeds the factor ^ 2 
difference expected from i = 70° . This is also discussed in 
subsection 5.5. 

Over the 30-200 keV range, the ratio in fig- 
ure 12 increases with energy the E roughly as cx 
GRO J1655— 40 has a harder continuum. This is con- 
sistent with the PL approximation to their hard X-ray 
spectra, which gave F = 1.75 for Cyg X-1 (section 3.1) 
and F = 1.39 for GRO J1655-40 (Paper I). In terms 
of the Comptonization scenario, this slope difference can 
be attributed to the higher y-parameter of compPSh in 
GRO J1655-40 (1.31+°;J^) than in Cyg X-1 (1.15^°;°^). 
This in turn can be explained rather securely by the higher 
Tc of GRO J1655-40, unlike the case of the HP vs. LP 
comparison of Cyg X-1 in which the origin of the change 
in y remained ambiguous. As discussed in Paper I and 
subsection 5.2, the higher To of GRO J1655— 40 may be 
due to its lower luminosity. 

In the PIN energy range, the ratio of figure 12 exhibits a 
slightly concave curvature. Because the continuta them- 
selves in this energy range should be simple power-laws 
dominated by compPSh, this effect must be caused by dif- 
ferences in the refiection features. In fact, figure 5 reveals 
two effects. One is that the reflection hump of Cyg X- 
1, relative to the Compton continuum, is inferred to be 
stronger than that in GRO J1655— 40. This is again con- 
sidered an inclination effect, because the reflected com- 
ponent is expected to be proportional to cosz while the 
primary continuum is thought to be rather isotropic. 
Actually, the values of U/27r in table 1, which is calculated 
after correction for cosi, becomes comparable between the 
two objects. The other effect is that the reflection hump 
in Cyg X-1, modeled with i = 45°, emerges at 30-50 keV, 
whereas that of GRO J1655-40, modeled with i = 70°, 
is peaked at higher energies. This is because a larger i 
makes the reflection more forward directed, and hence 
the reflected photons suffer less Compton degradation 
(Magdziarz & Zdziarski 1995). These two effects in com- 
bination are considered to produce the slightly concave- 
shaped spectral ratio in the PIN band. Incidentally, the 
comparable values of ri/27r between the two objects dis- 
favors the idea that the hot Comptonizing corona may 
have a semi-relativistic outflow velocity away from the 
disk plane (Malzac et al. 2001), because in that case we 
would observe a smaller value of n/2Tr from Cyg X-1 due 
to the enhanced continuum. 



To summarize, the spectrum of Cyg X-1 and 
GRO J1655— 40 arc both dominated by the Comptonized 
continua, but Cyg X-1 exhibits more prominently the flat- 
disk-rclated spectral features presumably due to its lower 
inclination. Another inclination effect manifests itself in 
the subtle differences in the 10-70 keV continuum shape, 
presumably due to inclination-dependent changes in the 
strength and peak energy of the reflection hump. These 
results arc consistent with the view that the Comptonizing 
hot cloud is geometrically thick, and hence its emission 
depends much less on the inclination than the flat-disk 
related components. 

Analyzing broad-band BeppoSAX spectra of 
XTE J1118-f480, Frontera et al. (2001b) discovered 
that this transient BHB has a very weak soft excess and 
a weak reflection hump. Like the present discussion on 
GRO J1655-40, these properties of XTE J1118-F480 
may also be ascribed to the inclination effects, since this 
object is also reported to have a rather high inclination 
as i = (68 ± 2)° (Gelino et al. 2006). 

5.4- The disk parameters in the low/hard state 

In the high/soft state of Cyg X-1, a diskbb fit to 
the ASCA data showed the optically-thick disk to have 
Hn = 71 km (Dotani et al. 1997), which translate via equa- 
tion (1) to i?in = 90 km for i = 30° (Makishima ct al. 
2000) or i?i„ = 100 km for i = 45°. The disk is securely 
considered to reach the last stable orbit (Dotani et al. 
1997; Makishima et al. 2000). In the LHS, a similar cool 
disk is considered to be present as well, but its conditions 
become significantly more difficult to estimate, because 
its emission is overwhelmed and strongly modified by the 
Compton process. As a result, there is considerable con- 
troversy over the innermost disk radius in the LHS. 

Theoretically, the Advection Dominanted Accretion 
Flows (Narayan & Yi 1995) and other related hot flow 
models, meant to explain the LHS, all assume that the 
cool disk is truncated at some radius larger than the last 
stable orbit. (One important difference is that the orig- 
inal ADAF model assumes synchrotron photons as the 
soft seed photons to be Compton up-scatterd into hard X- 
rays, whereas the present scenario assumes the cool disk 
to supply the seed soft photons.) A gradual inward shift 
in this truncation radius will give rise to the correlated 
spectral softening, increased reflection, and increased vari- 
ation power frequencies, as observed (e.g., Gilfanov et al. 
1999; Zdziarski et al. 1999; Zdziarski et al. 2004; Zdziarski 
& Gierhhski 2004; Done et al. 2007). When the disk 
reaches the last stable orbit, a transition to the high/soft 
state will take place. However, this view has been chal- 
lenged by some LHS observations which suggest the disk 
to extend down to the last stable orbit (Miller et al. 
2006b; Miller et al. 2006a; Rykoff et al. 2007). Here we use 
the broad-band Suzaku data to infer the disk parameters 
in the LHS. 

In subsection 3.1, the conventional diskbb+cutof f pi 
flt to the present LHS data tentatively gave r-^ ~ 13 km, 
and hence Rm = 18 km for i = 45°. This apparently implies 
that the cool disk extends down to a much smaller radius 
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than in the high/soft state. However, as pointed out be- 
fore (Ebisawa et al. 1996; Di Salvo et al. 2001; Ibragimov 
et al. 2005) and revealed by the present data (figure 5), 
the LHS spectra of Cyg X-1 and other BHBs are signifi- 
cantly more complex than is described by a single power- 
law. In fact, our final solution invoking double compPS has 
given an order of magnitude larger value, Rin = 250 km 
(table 1), which is consistent with the disk inner radius 
receding by a factor 2.5 compared to the high/soft state. 
The large difference in Rin between the two modelings for 
the same data can be explained in the following manner. 
In our final model, the flux in the 1-2 keV range is mostly 
explained by compPSs (figure 5a), with diskbb carrying 
softer residuals. If, instead, a single straight continuum 
is employed, the diskbb component becomes unusually 
hot trying to explain the 1-2 keV excess (Di Salvo et al. 
2001; Ibragimov et al. 2005), and hence the inner disk ra- 
dius becomes apparently too small. This gives a lesson 
that disentangling the disk parameters, where the disk is 
not the dominant component in the spectrum, is subject 
to large uncertainties in the underlying continuum shape. 
Evidently, the case of GRO 1655—40 is subject to larger 
difficulties because of the weaker soft excess. 

The Comptonization scenario involves another impor- 
tant effect which artificially reduces the disk radius: the 
disk photons are "lost" from the diskbb component into 
the dominant Comptonized continua. In estimating the 
true disk radius, to be denoted we must take this 

effect into account, knowing that Comptonization con- 
serves the photon number (Kubota et al. 2001; Kubota 
& Makishima 2004; Kubota & Done 2004). Practically, 
can be derived by quadratically summing R^i of 
diskbb, Rfl"'^ of compPSh, and of compPSs as 

{Rl^y^Rl + {R^^' + {R!^' . (5) 

Then, equations (3), (4), and (5) readily yield the total 
luminosity £disk of the cool disk as 

^disk = ^raw + -^scod = 47r (-Rf°*) CTo^lt • (6) 

From the values in table 1 for the time-averaged spectrum, 
we obtain ^ 330 km, which is equivalent to ^ 15i?g 
for Mbh = 15 Mq. The implications are that roughly half 
the disk area is directly visible, and the disk inner radius 
is factor 3 larger than in the high/soft state. 

Although the estimated value of may still be rela- 
tively small compared with the original idea of disk trun- 
cation, it can be subject to yet additional uncertainties 
which take part in when transforming the raw diskbb 
normalization into a physical radius. The stress-free in- 
ner boundary condition is probably appropriate for a thin 
disk extending down to the last stable orbit. In contrast, 
a truncated disk may have continuously rising stress to 
its inner radius. Similarly, the environment of the disk 
is different between the two states: in the LHS, the disk 
is likely to be more strongly heated by irradiation and 
conduction from the hot cloud. Therefore, the color hard- 
ening factor in the LHS may well be significantly larger 
than the value of 1.7 employed in equation (1); then, the 



true effective temperature would be lower, and hence i?in 
would become still larger. 

Three independent pieces of evidence further argue 
against the small disk radius in the LHS. One is the rela- 
tively small values of reflection, r2/27r~0.4, which requires 
the cool seed disk not to intrude too deeply into the hot 
corona. To reconcile the small values of n/27r with a con- 
figuration wherein the seed disk is nearly completely sur- 
rounded by the corona, we would have to invoke, e.g., mild 
relativistic coronal outflows (e.g., Malzac et al. 2001), but 
we already argued against it (subsection 5.3). Second, 
the intrinsic width of the Fe-K line, when interpreted 
as relativistic broadening in the disk near the BH, can 
constrain the innermost locations of the cool disk. For 
this purpose, we re-fitted the average spectra with the fi- 
nal model, but with the Gaussian replaced by so-called 
diskline model (Fabian et al. 1989). When fixing the in- 
clination at 45° , and assuming the line emissivity to scale 
as , the fit constrained the rest-frame line center en- 
ergy as 6.31q^ keV, and the innermost disk radius for the 
Fe-K line production as 13 jiy Rg which is fully consistent 
with the value of derived above. Finally, as quoted 
in subsection 4.5, Revnivtsev et al. (1999) studied fast 
variations of Cyg X-1 in the Fe-K line energy region, and 
argued the cool reprocessor to be located farther than 
^ 100i?g for Mbh = 15 Mq. Although the imphed ra- 
dius is much larger than considered here, this should be 
taken as another evidence supporting the disk truncation. 
(Incidentally, the present results on the fast variability in 
the Fe-K lines and the Compton humps constrain their 
production sites to be closer than 1.3 x 10^ Rg.) 

Based on these examinations, we estimate that the cool- 
disk radius of Cyg X-1 during the present LHS observation 
is ~ 15i?g ~ 330 km nominally, but is subject to rather 
large uncertainties, and is likely to be considerably larger 
than this. 

All the considerations presented here will also apply to 
previous LHS studies, where the small inner disk radius 
derived from simple continuum fits was used to challenge 
the truncated disk model (Miller et al. 2006b; Miller et al. 
2006a; Rykoff et al. 2007). In short, the small disk radii 
claimed for some BHBs in the LHS would need a careful 
revision. 

5.5. A possible interpretation of the fast X-ray variation 

The remaining issue is how to interpret our results on 
the time variation, summarized in subsection 4.6 and sub- 
section 5.1. In a standard accretion disk which we be- 
lieve is approximately correct in the present case, the disk 
parameters are mutually related as Tin oc {R\°^)~^^^M^/'^ 
(e.g., Makishima et al. 2000), where M is the mass accre- 
tion rate. Since Tin has been confirmed to change little 
as Cyg X-1 varies, the variation must also keep and 
M approximately constant. (An alternative case, namely, 
positively correlated changes in and M, would be 
unlikely, because an increase in M would enhance ra- 
diative cooling; and hence reduce -Rf°'; Zdziarski et al. 
1999; Zdziarski & GierUiiski 2004; Done et al. 2007.) As 
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a result, the total disk luminosity Ldisk, defined by equa- 
tion (6), should also be kept rather constant. However, 
Lseed has bccn found to increase from LP to HP (subsec- 
tion 4.4). Therefore, Lraw must decrease in HP by 20- 
30%. As confirmed in subsections 4.4 and 4.5, the data 
permit (though do not require) such a decrease in Lraw 
from LP to HP. 

Based on the above arguments, we presume that Cyg 
X-1 gets brighter as a larger fraction of the disk photons 
are intercepted by the Comptonizing corona, while the 
underlying cool disk itself is kept essentially unchanged. 
This will cause an increase in Lseed, and the associated 
decrease in iraw- Although this inference does not nec- 
essarily specify how the corona is actually changing, we 
may consider one possible scenario below. 

Let us, for example, assume that the Compton cloud 
has a rather spherical shape, or a thick torus geometry, 
with its outer radius located at Rcc ^ 100 Rg ^ 2 x 10^ 
km where the free-fall time is ~ 0.1 s. (The actual in-fall 
time scale of the corona may well be longer, depending on 
the viscosity.) Based on our discussion made in subsec- 
tion 5.4, the cool disk is considered to partially protrude 
into the corona, but truncated at, say, ^50 Rg. The 
Comptonization will then take place both in the overlap- 
ping region (between and i?cc), and inside of 
the corona. The corona is assumed to be rather turbu- 
lent with strong density inhomogeneities, so that we need 
two optical depths, Ts ~ 0.4 and Ts ~ 1.5, to express its ef- 
fects. Furthermore, we expect to increase as the corona 
falls toward the BH, for the following two reasons. One 
is because, compared to the corona inside Rm, that in 
the overlapping region must be more strongly cooled by 
the disk photons, and less efficiently heated by protons 
due to lower densities. The other is because, at about 
Rcc, the disk surface layer must evaporate into the corona, 
so it will take a fraction of a second (calculated at Red 
Guilbert et al. 1982) for the evaporated electrons to be 
heated by protons via Coulomb interaction. Then, the 
observed Tg ~ 100 keV should be considered as a radial 
average over the corona. 

Under the configuration as assumed above, we further 
speculate that, for some unspecified reasons, the corona 
has a number of "holes" (at least in the overlapping re- 
gion between the disk and corona), through which the disk 
is directly visible (subsection 5.2). With these, the spec- 
tral composition of Cyg X-1 (and of GRO J1655— 40) can 
be explained at least semi-quantitatively. Furthermore, 
this scenario can explain why the directly visible disk and 
the Fe-K lines of GRO J1655— 40 are weaker than those of 
Cyg X-1 beyond what is expected by the inclination differ- 
ence, because the holes would act as a vertical collimator 
to direct these disk-related emission components toward 
pole-on directions. The disk reflection would also be af- 
fected, but thier intrinsically broad shape would remain 
rather intact. 

In order to explain the fast variability, we may assume 
that the opening fraction of the coronal holes vary with 
time, due to some fluctuations in the disk evaporation 
rate. As the hole area diminishes by an increased evapo- 



ration, we expect Lsoed to increase and iraw to decrease, 
as observed in the changes from LP to HP. Since the incre- 
mented portion of the corona must initially have a lower 
electron temperature due to the finite ion vs. electron 
coupling time, the HP spectrum is expected to be less 
Comptonized than in LP, again as observed. Electrons in 
that portion of the corona will be gradually heated up, on 
time scales of a fraction of a second, as they fall inwards 
across This will explain, at least qualitatively, the 

XIS to PIN phase lag which would be too large to be ex- 
plained by time-of-flight delays of photons in the Compton 
process. 

6. Summary 

Through the detailed analysis of the Suzaku data of Cyg 
X-1, we have obtained the following firm results. 

1. The 0.7-400 keV Suzaku spectra are reproduced by 
a sum of a cool disk emission, two Compton con- 
tinua both with rcfiection {il / 2tt ^ A) , and a mildly 
broad Fc-K line. 

2. The two Compton continua, with t 0.4 and ~ 1.5, 
can be described by a common temperature of ^ 
100 keV, although the data do not exclude them 
having different temperatures. 

3. The disk is characterized by Tin ~ 0.2 keV and Rin ^ 
330 km ( ^ 15 Rg). Roughly half the disk emission 
gets Comptonized, while the rest forms the spectral 
soft excess. 

4. When Cyg X-1 brightens up on a time scale of 1- 
200 s, a larger fraction of disk photons becomes 
Comptonized, although the disk parameters remain 
unchanged. The spectrum softens, as a result of a 
decrease in either To or r. 

5. A comparison between Cyg X-1 (with i ~ 45°) and 
GRO J1655-40 (with i - 70°) supports the view 
that the disk has a flat geometry and the Compton 
cloud is more spherical. 

Combining these results with some theoretical consid- 
erations, we propose the following speculative scenario 
on the LHS of Cyg X-1 (and possibly of other BHBs). 
That is, the Componizing corona, with a roughly spher- 
ical (or thick torus) geometry, extends up to, say, ~ 
100i?g, and the cool disk protrudes half way into it. The 
corona is highly inhomogeneous (as required by the dou- 
ble Compton r), and has openings (as required by the 
soft excess) to allow a direct view of part of the disk. The 
hard X-ray intensity increases when the opening area of 
the corona decreases and hence a larger number of disk 
photons become Comptonized. 
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Table 1. Simultaneous fits to the XIS2, HXD-PIN, and HXD-GSO spectra of Cyg X-1 using two thermal Comptonization 
components.* 



Component 


Parameter 


Average 


High Phase 


Low Phase 


Difference 


GRO J1655-40 ^ 


wabs 




"•"-0.3 


(6.6 fix) 


(6.6 fix) 


(6.6 fix) 


' -^-0.1 


diskbb 


Tin (keV) 




0.19 ±0.01 


0.18±0.01 


(0.19 fix) 


0.18±0.01 




Rin (km) * 

-'-'raw 


250t-" 

^ ^ -I- 1 '^ 

iota 


250+^" 
.11 
1 1 + 


10±1 


< 110 
< 2 


< 73 

< 1 


compPSh 


To (keV) 
optical depth 
y-parameter 


100 ±5 
1 ^q+0.07 

-•^■^3-0.03 

75± 10 


^•^^-0.09 

n^A-l-l ^ 


105±ili 

1 4+0-4 
J^-^-o.i 

-0.03 

r^r'-^-^ 

751}" 


9 (-1+0.5 

0.89 ±0.3 

45 ±5 


140 ±10 
1.14±0.08 

- -)_n 1 1 
^■'-*^-0.08 

26+^ 


compPSs^ 


optical depth 
y-parameter 
(km) II 


"•^^-0.03 

2001^3° 


0.3811;;^^ 
0-28Io;ii 

2101^1; 


0.321^:^^ 

"■^ ' -0.03 

210t^° 


0.4±0.2 

0-16lo;o4 
130lJ§ 


n oc+o.Oi! 
"•^'J-0.03 

^ r»n-l-0 07 
« , 4-7 


n ' ^1 

rcncction 


^ (erg cm s~^) 


T. — ,1(1 '1 

^■^-0.3 

<30 


(0.4 fix) 
(0 fix) 


(0.4 fix) 
(0 fix) 


(0.4 fix) 
(0 fix) 


0.5±0.1 

400 ±100 




^in^(^g) 


2500^^ 


(2500 fix) 


(2500 fix) 


(2500 fix) 


> 200 


Gaussian 


centroid (keV) 


6.3±0.1 




6 3+"-^ 

D.<3_o.2 


(6.3 fix) 


6.3i[ii 




eq. width (eV) 


290 ±90 


290 ±50 


330tl° 


400 ± 200 


iooti° 




sigma (keV) 


^■"-0.2 


"•^-0.1 


1.1±0.2 


(1.0 fix) 


0.7±0.1 


fit goodness 


xt {^) 


1.13 (349) 


1.08 (228) 


1.17 (193) 


1.01 (64) 


1.03 (1484) 


figure 




figure 4a, d 


figure 9a 


figure 9b 


figure 9c 


Paper I 



* Cyg X-1 and GRO J1655— 40 are assumed to be at distances of 2.5 kpc and 3.2 kpc, respectively. The fitting 

model is wabs* (diskbb+compPS+compPS+gauss) . 
t The results on GRO J1655-40 refer to Paper I. 

t The innermost disk radius derived via equation (1). assuming i ~ 45° for Cyg X-1 and i = 70° for GRO J1655— 40. 
^ Bolometric luminosity in 10'^^ erg s~^, obtained as iiaw = ^'^Rfn'^oT^^ with (Tq the Stefan-Boltzmann constant. 
'I Square root of the compPS normalization, converted to the inner radius of the seed photon disk via equation (2). 

Note this is not divided by the vcosi factor. 

Assumed to have the same Tc as compPSh. 
% The parameters are assumed to be common to the two compPS components. The reflector inclination is assumed 

in the same way as in -I-. 
L Not constrained in the fitting. 



Table 2. Cutoff power-law fits to the HXD (PIN and GSO) spectra in HP and LP.* 



Component 


Parameter 


Free 


Common F 


Common i?cut 


cutoff PL 


r (HP) 


1.51 ±0.02 


1.51±0.01 


1.52±0.01 




r (LP) 


1.50 ±0.02 


_ t 


1.49 ±0.01 




£^cut (HP) 


212 ±12 


208+^ 


220 ±10 




Scut (LP) 


234 ±16 


240 ±12 


_t 




Norm, ratio-f 


1.31 


1.27 


1.39 


fit goodness 




0.93 (168) 


0.93 (169) 


0.94 (169) 



* The HP and LP spectra were fitted simultaneously with the cutoff pl+pexrav model. The reflection strength 

was fixed at f2/27r = 0.16 
^ Constrained to be the same between the HP and LP spectra. 

The ratio of cutoff pi normalization (at 1 keV) between HP and LP. 



